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Abstract 

The one-loop effective energy density of a pure 577(2) Yang-Mills theory 
in the Savvidy background, at finite temperature and chemical potential is 
examined with emphasis on the unstable modes. After identifying the stable 
and unstable modes, the stable modes are treated in the quadratic approxi- 
mation. For the unstable modes, the full expansion including the cubic and 
the quartic terms in the fluctuations is used. The functional integrals for the 
unstable modes are evaluated and added to the results for the stable modes. 
The resulting energy density is found to be real, coinciding with the real part 
of the energy density in the quadratic approximation of earlier study. There 
is now no imaginary part. Numerical results are presented for the variation 
of the energy density with temperature for various choices of the chemical 
potential. 



Keywords: Savvidy vacuum; 577(2) Yang-Mills; background chromomag- 
netic; finite temperature; chemical potential; unstable modes; quartic terms 
in fluctuations. 



1 e-mail : sar athy @imsc . res . in 
2 e-mail:alok@imsc.res.in 



1 



I.INTRODUCTION 



The ground state of SU(2) Yang-Mills theory in a covariantly constant 
chromomagnetic field in the third color direction (Ff 2 = H) had been studied 
by Savvidy [1]. The important result obtained was the one-loop effective 
energy density was lower than the perturbative vacuum and the theory has 
a gluon condensate, < \F^ v F^ va \ >^ 0. However, Nielson and Olesen [2] 
pointed out that the one-loop effective energy density had an imaginary part 
due to the lowest Landau level. Various studies were devoted to this aspect 
[3]. In recent times, gauge invariance arguments to exclude the imaginary 
part [4] and a stable vacuum by dynamically generating mass term for the off- 
diagonal gluons [5] have been proposed. In all these studies, only the terms 
quadratic in the fluctuations are kept, throwing out the cubic and quartic 
terms in the fluctuations. As early as 1983, the Savvidy vacuum within 
the framework of the background field method was examined by Flory [6] 
and later by Kay [7] treating the unstable modes including the cubic and 
quartic terms and the possibility of obtaining real effective energy density 
was pointed out. The present authors with Kay [8] have explicitly shown 
at zero temperature, the resulting energy density, when the cubic and the 
quartic terms in the fluctuations are included for the unstable modes, has no 
imaginary part and coincided with the real part of the earlier calculations. 
Briefly, the stable modes are treated in the quadratic approximation and for 
the unstable modes, the functional integral is evaluated keeping the cubic 
and quartic terms along with the quadratic term. This result (Eqn.30 of [8]) 
has no imaginary part and when added to the result of the stable modes, 
gives the effective energy density which is real. The calculations leading to 
this result are gauge invariant [8]. 

The Savvidy vacuum at finite temperature has been studied by many 
authors [9 - 15]. In all these studies, the quadratic approximation has been 
used. As a result, the effective energy density involved an imaginary part, 
dependent on the temperature. This is serious as it persists at high tem- 
peratures. It has been observed by Meisinger and Ogilvie [15] that with the 
introduction of a non-trivial Polyakov loop, specified by 0, it is possible to 
stabilize the vacuum, if (3yJgH < <p < 2n — P^gH, where /3 = 1/kT, H is the 
chromomagnetic background in the third color direction, as for this range of 
0, the imaginary part becomes zero. However, the imaginary part is non-zero 
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at the global minimum. So, in the understanding of the finite temperature 
behaviour, the imaginary part severely inhibits the progress. 

It is the purpose of this work to extend our method [8] at zero tem- 
perature, to finite temperature with chemical potential, by separating the 
unstable modes and including the cubic and quartic terms in these modes in 
the evaluation of the functional integral and adding to the contribution of the 
stable modes in the quadratic approximation. The motivation is to get real 
effective energy density. Another purpose is to resolve a discrepancy in the 
analytical expressions for the one-loop energy density between [14] and [15]. 
The discrepancy is an interchange of J\ and Y\ and also a relative sign be- 
tween the two Ki functions in the energy density. Our results show that the 
finite part of the effective energy density is real. There is no imaginary part. 
The real energy density coincides with the real part of [15] which therefore 
resolves the discrepancy alluded above in favour of [15]. 

The next section develops the formalism for the evaluation of the effective 
energy density. The stable and the unstable modes are separately treated. 
Section. Ill provides the details of handling the unstable modes including the 
cubic and quartic terms in the expansion. The results are added to the 
contribution of the stable modes and the full expression for the effective 
energy density is exhibited. Section IV contains the details of the numerical 
evaluation of the effective enerdy density for various temperatures. 

II. EFFECTIVE ENERGY DENSITY IN THE BACKGROUND 
FIELD METHOD (BFM) 

The Euclidean functional integral for an SU(2) pure YM theory is 



Z 




(1) 



where 



S 




(2) 



and 



F; v = d,A a u -d u A; + ge abc AlA:. 



(3) 
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Expanding A ^ = A ^ + with A® as the classical background field satisfying 
the equation of motion D ab F b v = 0, with D ab = d^S ab + ge acb A c ^ as the 
background covariant derivative, F® u is same as (3) with A® and using the 
background gauge 

D a X = °> ( 4 ) 

we have 

Z = j[da«y, (5) 

with 

- j{(«X) 2 -°X a X}-^^(-^?))' ( 6 ) 

where 

= (D?D b x c )5 lxu + 2ge aec F* u . (7) 

In arriving at (5) and (6), we have introduced the gauge fixing and the 
ghost lagrangian for the background gauge (4) and integrated the ghost fields, 
resulting in the last term in (6). The expansion in (6) is exact. If the cubic 
and the quadric terms in are neglected in (6), then the effective potential 
will be given by 

T 1 -^ = ^Treogdet(-e™)-Treogdet(-Dl b D b ;). (8) 



The Savvidy background is chromomagnetic in the third color direction 
and is given by 

A« = ; A1 = ^ 3 (-^,^,0), (9) 

which gives Ff 2 = H and solves the classical equation of motion, D ab F b u = 0. 
We introduce the chemical potential (/i) as a background field, following 
Loewe, Mendizabel and Rojas [16], namely 

K = -v,5 a3 , ^ = (1,0,0,0). (10) 
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This method, introduced in [16], has the advantage that as chemical poten- 
tials are not introduced as Lagrange multipliers in BFM, it is not necessary 
to compute the conserved charges. In [16], only the background (10) was 
considered. In our case we combine (9) and (10) as 

Ai = ,f ,o), (ii) 

which gives again F 3 2 = H as the non-vanishing background field strength 
and which solves the classical equation of motion. In this method of intro- 
ducing the chemical potential (10,11), it will play the role as the Polyakov 
loop specified by a constant A% field in the third color direction as in [15] 
with the identification of = fifi. As in [8] (eqn.10), we notice that (7) gives 

©44 = ©33 = DfD\\ (12) 

where now 

Df = d x 5 ah + ge a3b Al + fie a3b v x , (13) 

where A\ is given by (9) and v\ in (10). The relation (12) gives the important 
result that the contributions of 644 and 633 to Y 1 ~ loop in (8), cancel the ghost 
contribution in (8). Using (11) in (7), it is easy to see that 

e S = e « = =e- = e- = o. (u) 



The remaining eigenvalues are from 0^ c for i,j = 1,2 only. The eigen- 
modes and the corresponding eigenvalues are: 

(a\ + ia\) +i(a\ + ia 2 2 ) : (k 4 + fi) 2 + kj + (2N + l)gH - 2gH, 

{a\ + ia\) -i{a\ + ia 2 2 ) : {k A + fi) 2 + kj + {2N + l)gH + 2gH, 

(a\ - ial) + i(a\ - ia 2 2 ) ■ (h ~ v) 2 + k\ + (2N + l)gH + 2gH, 

(a{ - ial) - i(al - ial) (h ~ + ^3 + (2iV + l)gH - 2gH, 

a 3 : k\ + k\ + k\ + kl; i = l,2, (15) 

where iV = 0,l,2,---is the harmonic oscillator (in the x — y plane) quantum 
number. 
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We now pass on to the finite temperature case by replacing [17] k^ by 
and /f^ dki by J2^=-oo- With this, the effective potential in the "quadratic 
approximation" becomes 

1 i r ^ - - r dh 

2 ^2*)^ tJ~™ 2tt 



A 2LV /3 
1 27m 

+ ^— {(-^- + /i) 2 + A;3 2 + (2iV+ l) 9 H + 2gH} 
+ ^-{(^- / .) 2 + ^ + (27V+ l) 9 H + 2gH} 
+ iogj- 2 {(^-^ + kl + (2N+l)gH-2gH}), (16) 

where the first term corresponds to the last eigenvalue in (15) and the re- 
maining correspond respectively to the first four eigenvalues in (15). The 
pre-factor ^ is the harmonic oscillator degeneracy and the pre-factor 2 in 
the first term accounts for the two modes in the last term in (15) for % — 1, 2. 
A is the dimensionful parameter to render the argument of the logarithms 
dimensionless. We suppress this hereafter. 

This expression (16) agrees with Eqn.8 of [15] with their <f> replaced by fj,(3 
and with Eqn.(2.16) of Ninomiya and Sakai [12] with /j, — 0. It can be seen 
from (15), that for k$ = and n = 0; N = 0, the eigenvalues first and the 
fourth become fi 2 —gH and therefore to avoid negative eigenvalues /i > \J gH. 
However, for n — 1, the fourth eigenvalue becomes — fi) 2 — gH (with N=0 
and A; 3 = 0) and if this is to remain positive, then \i < ^ — \J gH. With this, 



the first eigenvalue remains positive. So for \fg~H < /i < y — \fgH, it is 
possible to avoid the negative eigenvalues and hence the instability. However, 
the instability enters at the global minimum [15]. We will return to this in 
the next section. Now, we proceed to evaluate (16). In order to exhibit 
the contribution from the unstable modes, we consider the second logarithm 
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(without its prefactors) in (16) and write that as [8] 

00 00 dk* r°° 



- E Ef - H dt f 1 e^^^ +A \ (17) 

n=-oo N=0 J -°° 2?F ^° 



where A = 2NgH — giL After performing the k 3 integration and the sum 
over n, we find 



^ = -^7=E/ tftr ' 3(-^-,-^l' 



where 



00 

inrn 2 „2ninz 



83(2, t) = E • 

n=— 00 



Now, using the property of the # 3 -function [18], 
e 3 (z,ir) = 



'ir' r 
we get 

?2 



12 - ~^hJo dtt H ^^l U 



N=0' 

Including the prefactor in (16), the contribution to r 1_ ' oop from L2 is 

N=0' 



The unstable mode corresponds to iV = 0. Splitting the above sum over iV 
as for N = and N — 1, 2, • • • , 00, we find 

r ^ (i2 ) = _^*^ ( £g )( ^ + _J^), (U!) 
where the first term in ( • • • ) is the contribution from the unstable mode. 

The third term in (16) has no negative eigenvalue and writing this as in 
(17), the sum over TV = 0, 1, 2, • • • , 00 is performed to give 

r^ (i3) = -^ 2 f dtt -^^ )T g^. ,19) 
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The fourth term in (16) is the same as the third term except for // replaced 
by —[/,. So also the fifth term and the second term are same but for this 
replacement. As 9 5 (z,t) = 9 3 (—z,r) [18], we have 



r i-ioo P(L , _ 9H I- 2 tP 2 \ ( tgH e-W + x 



8tt 2 



(20) 



Now ^(ff , f^) can be rewritten using 



<H--t) = J2 e "" ' 

£=—oo 



as 

?2 



The exponentials in (20) can be written as C g^h(tgH) an< ^ so (^0) becomes 

-S/o M ^" 2 ^«(l + 2E^co S (^) e-4f). (21) 

The first term in (16) is the familiar term and its finite part is 
Thus the one-loop potential is 



45/3 4 ' 



7T 2 

Y^—ioop 



45/3 4 



where r = g-f/t. The zero temperature part is reproduced. In [8], this was 
evaluated treating the unstable modes carefully including the cubic and the 
quartic terms and the effective energy density has no imaginary part. It was 
found to be, 

£i- loop (T = 0) = ^ + ^ L (tog(<L r )--). (23) 



We now proceed to evaluate the finite temperature part of (22), 



7T 2 



45/? 4 

(gHf /■«) _ 2 cos/i2r cx 



47T 2 Jo 

7T 2 



/ dr r — ) cos(fi[3£) e ^ ? 

Jo sinhr 



45/? 4 

(##) 2 /■«>, _ 2 , r _ r 2e- r . ... 

(24) 

where we exhibit the unstable mode contribution explicitly by the first term 
in (•••)■ Expanding 



oo oo 

2nr „—2t i \ ^ — 2nr 



1 e n=0 n=l 



e 



and introducing 

Ji = / drr- 2 e T e ~, (25) 

Jo 

I 2 = r drr~ 2 e-^, (26) 
Jo 

h = 2 r dr t- 2 e ^ 2n+1 ^ e -T, (27) 
Jo 

the finite temperature part of (24) is written as 

7T 2 

Y^—loop 



45/? 4 

^!g cos (^) (/ 1 + / 2 + / 3 ) ! (28) 

in which I\ is from the N = unstable mode. The integrals in (27) are to be 
evaluated. 
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In l\ (25), we perform a Wick rotation to arrive at 
h 



-i 



i ~ / dt t e 
2?r« 



where ii/j 1 ^ is the Hankel function of the first kind [19]. 
I 2 and J 3 integrals are evaluated using [19] 



as 



K v {xz) = - / r"" 1 e-fC+V) rff, 
2 Jo 



(29) 



(30) 



(31) 



'3 = 8£ 



(2n + 1) 

so that the finite part of (28) becomes 

J2 



pi— loop^rp^ 



71 



45/? 4 
(9H) 



X 



(32) 



(33) 



n=l 



Using H[ 1] {(3£y/g~H) = Ji {(3£VgH) +iY 1 (/3ty/gH) , we find that our result 
(33) agrees with Meisinger and Ogilvie [15] who used a different method to 
evaluate (16). The interchange of J\ and Y\ and a relative sign between the 
two Ki functions in Starinets, Vshivtsev, and Zhukovskii [14] are incorrect. 
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From (28) and (33), it is seen that the unstable mode (N = 0) contributions 
are contained in 1\ and explicitly 

(34) 

The imaginary part above is remniscent of the zero temperature situation. 
In the later case, we [8] have treated the unstable modes including the cubic 
and quartic terms in (6) and showed then the contribution is real. In view of 
the important difficulties raising from the imaginary part at finite tempera- 
ture, we, in the next section, treat the unstable modes at finite temperature 
including the cubic and the quartic terms in the expansion (6). 

III. FINITE TEMPERATURE UNSTABLE MODES-INCLUSION 
OF CUBIC AND QUARTIC TERMS 

There are two unstable modes when the harmonic oscillator quantum 
number N = 0, as seen from (15). They are the first and the fourth line 
in (15) with N = 0. As the expression (16) was obtained in the "quadratic 
approximation", we cannot use (16) for the unstable modes when we want 
to include the cubic and the quartic terms in the unstable modes. So, we 
take the equation (5) for Z but confine ourselves to the unstable modes only 
here. Without loss of generality we will take \x = 0, zero chemical potential. 
The normalized unstable mode eigenfunctions are 

KM = ^§^ (xl+4) -JZtf e~^ + ^\ (35) 

where the first exponential is the ground statet wavefunction (JV = 0) of 
the harmonic oscillator in the (xi — x 2 ) plane, the second exponential is the 
box- normalized plane wave in the x 3 direction and the ^-harmonics in the 
^-direction. The index n is the Matsubaro index and is not the same n in 
(33) as the later index originated in the expansion of 1 _ e 1 _ 2T in the expression 
above (25). The unstable eigenmodes in (15) are then c(k 3 , n)(j)k 3t n(x) and 
c(k 3 , n)*(j)* k3jn (x) respectively. 
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The unstable modes involve Lorentz indices 1 and 2 (as can be seen in 
(15)) and the SU{2) indices 1 and 2, since the classical background in (11) is 
in the third isospin direction. Consequently, the cubic term e acd (D^ e aJ 1 ) a c a d v 
vanishes as in [8]. The terms quartic in in (5) are simplified to | a\ aZ a\ aZ 
Then, the full partition function for the unstable modes, from (5) and (6) is 

Z unstable = J [daft e -fd 4 *{au(q + kl- 9 H)a u+ £at}^ (gg) 

where a u stands for and aZ- From (36), the unstable modes for h\ + 
(^p) 2 < gH render the quadratic term in a u in the exponent of (36) 
divergent. However, the quartic term in a u in the exponent of (36) provides 
the necessary and the crucial convergence. Thus the overall integral over a u 
will be convergent. Now expanding a u in terms of the eigenfunctions (35) 
and carrying out the d 4 x integration in the exponent, the quadratic term 
becomes 

{k! + C-f) 2 -gH}c 2 (h,n), (37) 

and the quartic term becomes 

6^^(k 3 ,n) = f^c\h,n), (38) 

where we have taken all the four (f>k 3 , n ( x ) having the same fc 3 , n. Then, we 
obtain 

Z unstable = (/n fcB>B dc(Jbs,n) e ^-^ 3 ?>V-^}) D > (39) 

where D is the degeneracy factor D = with V as the spatial volume. 

Now, introducing c = \fgHc and k 3 = -Jl=^j, (39) is written as 

ZmtM , = ( i j^ja^-^-^y ,40) 



The contribution of the complete unstable modes to the energy density 
{-^v tog Z unstoh i e ) is then 

-^^JdhEn togJ(k 3 ,n), (41) 
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where 

dce i{ 3 W> ] >w> c} . (42) 

-oo 

In (42), if the quartic term is neglected, then (41) will give (34). We retain 
the crucial quartic term contribution. The integral over dc is convergent 
irrespective of the sign of the coefficient of the c 2 term. It is evaluated using 
[19] 



with 



/? 2 = ^ 



128tt(##)' 
Then (42) becomes 

x ^(«^Q { ^ + ( ^_ 1} .). (43) 

When fc| + { ^g H ) 2 < 1 (that is where the instability arose), the above 
expression using [20] (since n starts from 1, the argument of Ki will be small) 

K v (x) — » ; small x, 



becomes 



2-!r(i) ^ 



{*l + (^) 2 -i}* 2(7T^)i' 
^(vr^)32-fr(i). (44) 
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IS 



In particular the radical + (j^=[j) — l} 2 gets cancelled. This result i 
real. The imaginary part coming from the radical gets cancelled by the contri- 
bution from Ki . This is made possible by the inclusion of the quartic term. 

4 

When (44) is used in (41), the /^-integration and the sum over n (constrained 
by fc| + { p™g H ) 2 < 1) produces an uninteresting term which is omitted. 

When k\ + (jj==) 2 > 1, the expression (43) takes the form, using [20], 



/ 7T _ 

K u (x) -> Y^ e_:E ; large x, 



as 



J(k,n) = {kl + {^=?-l)\ (45) 
Then (41) is evaluated as: 

J dk^ogJikn) = -\rEt°9{% + (^^) 2 -l}, 

= ~\ J™ dk 3 iog{co S h(^(k 2 3 - l)gH) - 1}, 

= ~\ j~ 'dk 3 (M \kl - l)gH + 2eog(l-e-Py/®-V* H j), 

1 r°° - 7r 00 1 / — 

= —Q\ dk^{kl-l)gH - J2~ Y ^\/9H) + ^ 

(46) 



where 



h = 



gHY [ 2 cos{nf3\fgHcose)cos9 d9. (47) 

n=l J ° 



The first term in (46), when used in (41), produces a /^-independent 
contribution and neglected for finite temperature effects. The integral I\ in 
(47) is evaluated [19] to be 

T /^W f 1 i V- {-l) k ~\nP^g-Hf k x 



k=i 
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which is not contributing to the finite part of the energy density. Thus, the 
unstable mode contributions to the finite /9-dependent part of the energy 
density is found to be 



gH^/gH w v-^oo j_ 
2tt 2 /3 ' 2 ^n=l n 



(48) 



which is just the real part of (34). There is no imaginary part. This is due 
to the inclusion of the cubic and the quartic terms in the unstable modes. 
Thus, the difficulties associated with the imaginary part are not due to the 
intrinsic property of the SU (2) chromomagnetic ground state but due to the 
use of the gaussian approximation. This reaffirms our earlier study [8] of the 
same system at zero temperature. 

The complete expression of the energy density of the SU(2) chromomag- 
netic state, including the zero contribution from [8] is 



The finite temperature part agrees with (the real part of) [15]. A similar 
result was found for the zero temperature case in [8]. 

IV. NUMERICAL RESULTS AND DISCUSSION 

The expression for the effective energy density (49) involves summation 
over £ and in most studies, the high and low temperature behaviours were 
examined. Even then, the summation is very involved as can be seen from 
[14] and [15]. As the high and low temperature expressions are given in [14] 
and [15] (also in 12]), we attempt to evaluate (49) numerically. First, it can 
be seen from [21], the K\ functions have a fall-off to zero when the argument 



£ 




+ 




n=l 
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is greater than 5. On the other hand the Y\ function is oscillatory with 
decreasing amplitude. Second, we used the "polynomial approximations" for 
these functions given in [21] and verified that they are good by computing 
these functions for various values of x from 0.05 to large values and comparing 
them with the tables of these functions. Then, we used MATLAB to find 
the values of Yi(x) and Ki(x) for various values of x which agree with the 
previous method. Third, we need to evaluate the sums in (49). 



For the Y 1 function appearing in (49), we found the typical sum 



converged to a steady value for £x upto 200 whereas for the 
steady value is reached when £x = 20. These allow us to choose £ max = 200/a; 
for the Yi sum while i max = 20/x for the sum involving K\. Keeping £ max 
as 20/x, in the last K\ sum, the n sum was carried out from n — 1 to n max 



with n max = ■^{- ! fr i — 1}. In order to evaluate the temperature variation of 
(49), we first set (3 = and /j, = b^JgU. Then the temperature dependent 
part of (49) becomes 

£ T 7T 2 1 ~ COs(aM) / 7T 



(gH) 2 Aba 4 n 2 af^ £ ^ 2 

oo 

+ K x {at) + 2 V2n+lK 1 (aeV2n+l)y (50) 



In Fig.l, we have plotted j^js with T = ^, that is, T in units 

of ^W-, for b = 0,1,2,3. For 6 = 0, zero chemical potential, the variation 
is smooth apart from small oscillatory behaviour at low temperatures. For 
6 = 2,3, the variation shows a minimum and then raising smoothly. At high 
temperatures, the behaviour is like that of non-interacting relativistic gas. In 
[15], the Polyakov loop is measured in terms of which in our notation is /ifl 
and that is ab. In the sense that a is varied, it is not possible to relate directly 
our results to [15]. However, the importance of the chemical potential is seen 
in Fig.l. A non-zero chemical potential or non-zero <p triggers a possible 
deconfmement phase transition. Our variation is qualitatively in agreement 
with [15] for their "real part". 

Now, we wish to examine the inclusion of the cubic and quartic terms 
for all the modes. It can be seen from (15), the stable eigenvalues being 
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distinctly different from the unstable eigenvalues for a given N. So, it is 
justifiable to consider the corresponding eigenmodes as orthogonal. Then, 
from (15), it follows that the cubic terms will vanish for the stable modes as 
well. The resulting full expression can be evaluated as in (43) with explicit 
N appearing. When the logarithm is taken, as in (41), the finite part will 
remain unaltered. The situation for the unstable modes is different in the 
sense the troublesome imaginary part does not appear in (44). 

To summarize, we have considered the one-loop effective energy density 
of a pure 577(2) Yang-Mills theory in the Savvidy background, at finite tem- 
perature and chemical potential. The unstable modes are treated by keeping 
the cubic and the quartic terms in the fluctuations. This result is added to 
the contribution from the stable modes. There is no imaginary part. The 
variation of the energy density for a given chromomagnetic background with 
temperature is studied numerically. When the chemical potential is non-zero, 
the variation shows a minimum which is (roughly) interpretted as indicating 
a deconfinement phase transition. At high temperatures, the behaviour is 
like that of a relativistic gas. 
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Figure 1: Variation of scaled energy density with scaled temperature 
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